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The entanglement spectrum, i.e., the full distribution of Schmidt eigenvalues of the reduced 
density matrix, contains more information than the conventional entanglement entropy and has been 
studied recently in several many-particle systems. We compute the disorder-averaged entanglement 
spectrum, in the form of the disorder-averaged moments Trp^ of the reduced density matrix pa, for 
a contiguous block of many spins at the random-singlet quantum critical point in one dimension. 
The result compares well in the scaling limit with numerical studies on the random XX model and 
is also expected to describe the (interacting) random Heisenberg model. Our numerical studies on 
the XX case reveal that the dependence of the entanglement entropy and spectrum on the geometry 
of the Hilbert space partition is quite different than for conformally invariant critical points. 

PACS numbers: 



I. INTRODUCTION 

The ground state of a system at a quantum critical 
point shows universal behavior in many quantities. Cor- 
relation functions, for example, show universal power-law 
behavior, and in some cases these power laws can be ob- 
tained exactly by mapping the quantum critical point 
to a system in one more dimension. The most power- 
ful example of this mapping is for one-dimensional (ID) 
quantum critical points (QCPs) that become 2D classi- 
cal critical points with conformal invariance. In addition 
to standard correlation functions, it is now understood 
that the entanglement entropy, reviewed below (see the 
comprehensive reviews in Ref. [T]), is universal at such 
quantum critical points, and determined by the central 
charge of the associated 2D conformal field theor y^HU p or 
a partition of an infinite ID system into a finite chain of 
length I and the remainder, the entanglement entropy 
(the von Neumann entropy of the reduced density ma- 
trix pa) for £ much larger than the short-distance cutoff 
a is asymptotically 

c I 

Svn = -Tr [p A In p A \ = ~ In - + c, , (1) 
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where c is the central charge and c'j a non-universal ad- 
ditive constant. 

Other properties related to entanglement are less well 
understood, even at these quantum critical points, such 
as the entanglement spectrum (the full set of reduced 
density matrix eigenvalues) and the full set of entangle- 
ment Renyi entropies; one exception is free Fermi mod- 
els, where the entanglement spectrum is given by the 
spectrum of an effective "entanglement Hamiltonian'^. A 
form for the spectrurrP at ID conformal QCPs that is ex- 
act in some cases and a good approximation in others^' 
can be used to develop a theory of how finite entangle- 
ment perturbs criticality in numerical studies^lAAJ The 
entanglement spectrum has also been applied to under- 
standing gapped (non-critical) topological phases^HS 



where it contains information about the edge excita- 
tion spectrum that goes beyond the universal constant 
in the entanglement entropjffMm The same is true for 
quantum 2D models with conformal invariant ground- 
state wave-functions^. Also results for a critical non- 
conformal ID model are available^. 

This paper studies the entanglement spectrum at 
"random-singlet" ID QCPs, in which quenched disorder 
leads to an RG flow to infinite randomness. We obtain 
the disorder-averaged moments of the Schmidt eigenvalue 
distribution analytically and compare them to numerical 
results on a special case with a free-fermion representa- 
tion, the random XX model. While these critical points 
are not conformally invariant (after mapping to a 2D 
problem, the imaginary-time direction has no random- 
ness and is hence very different from the spatial direc- 
tion), their disorder-averaged correlation functions have 
nevertheless been understood in many cases^HUl by real- 
space renormalization group methocP^. The entangle- 
ment entropy at random-singlet critical points was al- 
ready knowrp3H2U to show universal behavior similar to 
that at ID conformal QCPs [Eq. Q], with a modified 
prefactor of the logarithm (analogous to c) that was ini- 
tially viewed as an effective central charge for random 
systems. 

However, the results presented here indicate that this 
similarity does not extend to the full entanglement spec- 
tra, which are rather different. We start by considering 
the disorder-averaged Renyi entropies 



S a =— lnTr^], (2) 

where the bar denotes the average over quenched disor- 
der. These Renyi entropies S a are quite simple in the 
random-singlet phase: they depend only on the mean 
number of singlets across the partition used to define 
the entanglement, just as does the entanglement entropy. 
The Renyi entropies already behave differently than in 
the conformal case. However, in disordered systems S a 
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is not the right quantity that determines the entangle- 
ment spectrum via Laplace transform in cPl To obtain 
the averaged moments of the distribution, one should in- 
stead consider the entropies corresponding to averaging 
the disorder before taking the logarithm 
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This definition has also the advantage to maintain the 
relationship of the pure system between the Tsallis 29 
entropies (Tr [p^] — 1)/(1 — a) and the Renyi entropies. 
These moments of the entanglement eigenvalue distribu- 
tion reveal the full distribution of the number of singlets 
crossing a boundary and require an improved calculation. 
Both generalized entropies reduce to the Von Neumann 
one for a — > 1 
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The entropies S a and S a together with other proper- 
ties are then studied for the random XX model and the 
validity of our results is discussed for general random- 
singlet ground-states. The manuscript is organized as 
follows. In Sec. [IT] we present the random-singlet pic- 
ture and we derive the entropies S a within strong disor- 
der renormalization group. In Sec. |III| we introduce the 
probability distribution of singlet formation and use it to 
derive the entropies S a . Numerical tests of the predicted 
entropies and the discussion of their universality are de- 
scribed in Sec. |IV| Finally in Sec. [V] we report our main 
conclusions. 



II. RANDOM-SINGLET PICTURE OF THE 
RENYI ENTROPIES 

The ground state of a strongly disordered s = 1/2 
Heisenberg chain or of the disordered XX chain 
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is described by the random-singlet phase (RSP) for essen- 
tially any probability distribution P(J) of the coupling. 
When a system reaches this phase the ground state be- 
comes almost factorized in singlets between spins at ar- 
bitrary large distances. The configuration of the singlets 
depends on the coupling constants Ji , but several univer- 
sal properties emerge in the average over disorder that are 
independent of the disorder distribution itself. The phys- 
ical properties of a system in the RSP can be attained 
in an indirect way, i.e. without referring (manifestly) 
to the particular Hamiltonian. The real-space renormal- 
ization group approach (RSRG) is based on the picture 
that the strongest bond gives rise to a singlet, and the 
near-neighborhood spins can be described by means of 
an effective interaction from second-order perturbation 
theory. 



Considering the XX Hamiltonian ([5|, the Ma- 

Dasgupta rule^l for the effective coupling constant after 
a decimation, i.e. the formation of a singlet, is 
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where Jm is the strongest bond of the chain of size L 
and J i {J r ) is the near-neighborhood left (right) coupling 
constant. One of the most important consequences of ^ 
is that the distribution of the couplings after a sufficiently 
large number of decimations m, with 
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is substantially independent of the initial distribution: 
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where T is the RG flow parameter = In ■ J ^. The 

distribution ([8]) is the key to physical characteristics of 
the random-singlet phase. It is also the main ingredient 
for investigating the entanglement of spin blocks. In fact, 
for a spin block of length £ in a given RSP configuration 
with n singlets linking the spins inside the subsystem 
with the spins outside (which we call in-out singlets) the 
reduced density matrix is 
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Thus, the entanglement of a subsystem of size £ depends 
only on the mean number n of in-out singlets. In par- 
ticular the entanglement entropy, as well as any Renyi 
entropy ([2JI, is proportional to the number of in-out sin- 
glets 
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(This result has been also discussed in Refs. [2"8"landl3lTl) 
Ref . I^Hshows that the averaged number of in-out singlets 
can be deduced directly from the flow equation for the 
distribution of couplings fii 
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After some manipulation, this equation leads tcP^ 
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and so the entanglement entropy of a block of length £ is 



(13) 
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that calls to mind the behavior where to keep the notation compact we defined 



with a weight-factor ., 
in the absence of disorder with an effective central charge 
In 2. 

Consideration of the Renyi entropy rather than the 
standard entanglement entropy suggests that the similar- 
ity between the entanglement entropy with and without 
disorder is only superficial. Indeed in the RSP all Renyi 
entropies scale in the same way ( 10 1 . If we wish to define 



an effective central charge, we could use any conformal 
Renyi entropy^ 
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as starting point, so that the effective central charge 
would have any value in the range [(In 4) , (In. 2) ] 
while a runs from 1 to infinity. Also the central charge 
of the clean system c = 1 belongs to this range, mak- 
ing questionable any attempt to generalize the Zamolod- 
chikov "c-theorem'ES. This picture from Renyi entro py 
is consistent with the previous counterexample^ 31 * 33 * in- 
dicating that there is no version of the c-theorem for en- 
tanglement entropy that would describe the flow from 
clean to random systems^ or within random systems^. 

The disorder-averaged Renyi entropies at random 
quantum critical points are universal and already indicate 
that the random-singlet phase's entanglement is quite dif- 
ferent from the universal entanglement at ID conformal 
QCPs. However, since they depend on the same quantity 
(mean number n of in-out singlets) as the entanglement 
entropy, they do not probe new features of the random- 
singlet picture. In the next section we consider additional 
quantities that are sensitive to new features and directly 
probe a memory effect in the RSRG flow, or "repulsion 
between decimations" in RG space, that was a key step in 
obtaining the correct value of n. Numerical tests of the 
predicted Renyi entropies are described in Section |IV| 



III. GENERALIZED ENTROPY AND THE 
PROBABILITY DISTRIBUTION OF SINGLET 
FORMATION 

The disorder-averaged Renyi entropy in the RSP only 
reflects the averaged number of the in-out singlets. Thus 
it is not a natural measure of the full in-out singlet distri- 
bution P(n). or the probability distribution of the Renyi 
entropy. P(n) can be examined by considering S a in Eq. 
|3]). In fact, denoting with g(t) the cumulant-generating 
function of the in-out singlet distribution P(n) 



g(t)=\n(e nt ) = lnJ2P(n)e nt , 
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it is straightforward that 
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t = t(a) = (1 - a) In 2. 
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Through all the paper t will always denote this quantity, 
even when the a dependence is not specified. S a does 
depend on a in the RSP, unlike the Renyi entropy S a . 
We require Eq. ( 15 1 to not blow up when n — ¥ 00, and so 



(assuming a reasonable P(n)) we need t < correspond- 
ing to a > 1. We not not discuss a possible analytic 
continuation to a < 1 (that also in some clean systems 
can be complicated^) . 

From the real-space renormalization group (RSRG) 
point of view, singlets form at a constant rate with re- 
spect to a "RG time" /x, and this rate determines the 
logarithmic scaling of entanglement entropy. En route to 
calculating this rate, Ref. [22] obtains the expression for 
the distribution of waiting times for a decimation across 
a bond since the last decimation: 



3+\/5 , 
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The above distribution has been deduced neglecting non- 
universal terms coming from the starting disorder distri- 
bution: Eq. (18 1 is only asymptotically true. For ex- 



ample, we expect that the additive constant of the von 
Neumann entropy Svn should be disorder dependent. 

During the RG time between two decimations several 
processes can happen. The most probable one is the for- 
mation of isolated singlets. Considering only this process 
leads to the renewal equation 

(e"%= / d/i'/OO+e* / d M 'iV)(e"V M < ■ (19) 
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This equation can be solved by Laplace transformation. 
Calling f(s) the Laplace transform of /(//) 
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and in particular n = lim g (t). 
After simple algebra, we obtain 
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that via Eq. ( 16 1 gives S a in the RSP. This is the main 



analytic result of this paper. It is useful to rewrite it in 
terms of the mean number of singlets as 



g(t) = tA t n + tB t . 



(23) 
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where the multiplicative t factor is introduced to write 
more compact formulas for S a via Eq. (16 1. The two 
constants A t and B t are obtained by plugging ( 23 1 into 
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Notice that in Eq. ( 23 1 all the dependence of g/^\ (t) on 
fi is encoded in n. Tn this way, we also separated the 
universal In £ behavior (we remind n oc In £) given by A t 
from the constant one B t . We will come back to the 



discussion of the universal features of Eqs. (|23| and (p4|) 
in the next section when comparing with t 
results. 
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IV. NUMERICAL RESULTS 

In this section we present numerical evidences confirm- 
ing the critical scaling of the quantities calculated ana- 
lytically by means of RSRG. We also present results for 
which we do not have yet any theoretical explanation, 
like the finite size scaling in the RSP. 

The entropies S a and S a can be directly calculated for 
the disordered XX chain ^ , by generalizing the method 
of Laflorencie^l In fact, for any realization of the dis- 
order (i.e. any distribution of the bonds J;), the XX 
model can be mapped into a free-fermionic Hamiltonian 
by the Jordan- Wigner transformation cj = Ylj<i a j a 'i'' 
that leaves the eigenvalues of the reduced density matrix 
of a single block unchanged, because the transformation 
is local inside the block. Defining the correlation matrix 
C/n = { c \ c n)i the reduced density matrix of a spin block, 
that goes from the site lp + 1 to Iq + £, is the exponential 
of a free-fermion operatoJS12SI and it is completely charac- 
terized by the I x I correlation matrix C in which indexes 
run from Iq + 1 to Iq + £, that we call cf °' . The entan- 
glement entropy of the block in this configuration of the 
disorder is then given by 
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while the Renyi entropy is 
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where we stressed the dependence on the disorder config- 
uration ({J;}) and on the first site of the block Z + 1. In- 
deed, on a single realization of the disorder, translational 
invariance is explicitly broken. Only after taking the 
disorder average translation symmetry can be restored. 
Having the Renyi entropies for a single realization, al- 
lows to obtain the asymptotic results for the disordered 



model, by averaging over a large enough number of con- 
figurations (generated randomly according to the specific 
rules for {</;}). S a and S a are obtained by averaging S a 
or e( 1 ~ a * >Sa respectively. 

The method we presented is an ab-initio calculation of 
the Renyi entropies for disordered spin chains valid ev- 
ery time the model has a free-fermionic representation 
(as in XX or Ising chains). It is however numerically de- 
manding. A more effective numerical technique exploits 
the RSP structure of the ground-state. Starting from a 
given disorder realization, we construct a singlet where 
the strong bond lies and we proceed to decimation ac- 
cording to the rule in Eq. We repeat this procedure 
until we spanned all the chain. At this point we are left 
with a collection of singlets, and then, counting number 
of singlets connecting the inside of the block with the 
outside, we have the configurational Renyi entropies from 
the relation Sl o] ({Ji}) = ({ J/}) In 2. As for the ab- 
initio calculation, S a and S a are obtained by averaging 
over the disorder. Note that S^ SP does not depend on 
a by definition, since for any configuration S a = nln2. 
Oppositely S a depends on a because the average is taken 
over e ( 1_Q ) s ° and indeed some results for S a have been 
already reported^ by using this method. For complete- 
ness, we give few general features for an intuitive picture 
of the entanglement in the RSP. After a decimation (JfjJ, 
the renormalized bond is strongly suppressed, i.e. sin- 
glets repel. The singlets that stay inside the block in- 
volve always an even number of spins, thus the parity of 
the block gives the parity of the number of in-out sin- 
glets. The spins belonging to the longest bonds crossing 
the two ends of the chain can be also thought as bound- 
aries of two open chains. This suggests that in the RSP 
(as it is the case for clean systems^) the entanglement 
entropy of a block of I spins in a periodic chain is equiva- 
lent to twice the entanglement entropy of £/2 spins in an 
open chain with the block starting from the boundary, 
i.e. S'P criodic (^) « 2S'° pen (^/2). However, this argument 
does not provide information about the additive constant 
(in clean models, the difference of the two const ant ter ms 
gives the Affleck and Ludwig boundary entropy ^ * 37 * 38 *). 

To avoid confusion between the two determinations of 
the entanglement, in the following we will always refer to 
the first method as ab-initio while to the second as RSP. 
We stress that the RSP technique can be applied to any 
model with an RSP ground-state, as for example the dis- 
ordered Heisenberg chains or spin-1 chains^, while the 
ab-initio one only to models having a free-femionic rep- 
resentation. However, the ab-initio method has the ad- 
vantage to be exact by definition. Instead, by counting 
the number of singlets, we make the assumption that the 
ground-state has an RSP structure and that all the uni- 
versal entanglement physics can be extracted from this. 
Although both assumptions sound reasonable, it is al- 
ways worthwhile to perform in parallel the two numeri- 
cal studies. In fact, the numerical counting of singlets is 
not the same as the analytic expressions derived in the 
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Figure 1: Ab-initio Renyi entropies for a disordered XX chain 
of 1024 spins. The average is over 73000 realizations. The 
variation of the color shows results from a = 1 (upper line) 
to a = 2.9 (bottom line). The yellow line is the asy mpt otic 
Von Neumann entropy (a — 1) obtained by Laflorencie^. 



previous sections because, in order to provide analytic 
results, few further assumptions have been made (e.g. 
considering only the formation of isolated singlets etc.). 
In case of disagreement between formulas and numerics, 
making the two computations in parallel helps to under- 
stand if the error is in the approximations made to solve 
the equations or in the RSP assumption itself. 

A possible generalization (that is currently under 
investigation's ) is to understand if the RSP structure 
catches the entanglement of two disconnected blocks. 
This can be achieved by calculating ab-initio Renyi en- 
tropies (indeed there are not known formulas for the en- 
tanglement entropy when the subsystem consists of more 
than one spin block, only some expressions have been re- 
cently found for the first integer Renyi entropies^) and 
comparing with the in-out singlets from both block s. It 
has been shown for conformal critical models^EH that 
the entanglement of two blocks provides much more in- 
formation about the conformal structure than the single 
block one, and it is then worth investigating this issue 
also for the random case. 



A. Analysis of S a 

We computed ab-initio the averaged Renyi entropies 
S a for many different system sizes. In Fig. [T] we re- 
port the result for a chain of L = 1024 spins for the 
disorder average over a sample of 73000 realizations. For 
1 <C I <C L, the various curves are parallel, with the 
slope predicted by Eq. (13), i.e. the leading term of S a 



is a independent. The non-universal additive 0(1) term 
clearly depends on a, as in the clean case. On top of 
a smooth behavior, we can see oscillating contributions, 
evident for small I and large a. Their presence does not 



Figure 2: Ab-initio Renyi entropies for a disordered XX chain 
of 1024 spins minus the RSP value. The averages are over the 
same sample of 73000 realizations. 



come unexpected: also in clean chains^"^, there are os- 
cillating terms that (in zero magnetic field) are parity 
dependent, i.e. they are of the form (— However, 
for random systems the oscillations have a different form 
and they decay rather quickly with £ (as opposite to S a 
as we shall see) . We do not have a proper theory for their 
origin, nor a phenomenological description, but their un- 
derstanding is beyond the goals of this paper since they 
do not influence the determination of the asymptotic ba- 
havior. When I approaches the chain length L, sizable 
finite-size corrections are visible. Next subsection will be 
devoted to their accurate study, while here we continue 
with the asymptotic analysis of S a . 

We compare the data in Fig. [TJfrom the ab-initio cal- 
culation, with the numerical results obtained using the 
RSP approach on the same random sample of 73000 re- 
alizations of Ji. According to Eq. (10 1, the RSP Renyi 

For this 



entropies do not depend on a by definition, 
reason, in Fig. [2] we report the difference between the 
RSP Renyi entropies and the ab-initio ones presented in 
Fig. [T| After a transient behavior for small £, all the 
curves with varying a approach a constant, indicating 
not only that the universal leading logarithmic term in 
S a is correctly described by RSP, but also the finite size 
corrections are. In the range of a considered in the fig- 
ure, we find that the additive constant is well described 

by 



S a ~S* SP + -+6 + 0(1), 



(27) 



where the disorder-dependent constants a and b in the 
case of random disorder take the values a 0.61 and 
b « -0.47. 
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Figure 3: The finite-size scaling function for the entanglement 
entropy Y(x) in Eq. (291. Main: RSP data averaged over 
1440000 disorder realizations for L — 1024. The continuous 
(red) curve is the proposed phenomenological formula (311 
describing perfectly the data points. Inset: The same plot for 
different values of L, showing the collapse on a single scaling 
function. 



B. Finite-size effects 

Having established the correctness of the asymptotic 
RSRG results for S a in the region 1 -C I -C L, we can 
consider the finite-size effects. One of the most remark- 
able result of conformal invariance is that the finite-size 
scaling is obtained with the replacement 



£ — > — sin — 
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in the thermodynamic limit result. The rhs of the above 
equation is known as chord length. However, when con- 
formal invariance is broken, the chord length does not 
give the finite-size scaling. In fact, using the results re- 
ported above, it is easy to show that this is the case, as 
it was already shown for some random Ising systems^. 

Even if conformal invariance is broken, scale invariance 
still holds. Thus the finite size scaling can always be 
taken into account by the substitution 
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The great predictive power of conformal symmetry is that 
independently of the observable (but built with primary 
operators) the scaling function is always Y{x) = sin(x), 
while in general scale-invariant theories the function Y(x) 
does depend on the observable. Some results on the finite- 
size scaling of entanglement in ID critical non-conformal 
systems have been already reportecP^^. The function 
Y(x) for S a must however satisfy simple symmetry con- 
straints. First, S a is symmetric for £ — > L — £, thus 
Y(x) = Y(tt — x). Second, periodic boundary conditions 
require S a to be a periodic function of £ of period L, and 



so Y{x) = Y(tt+x). Thus we can expand Y(x) in Fourier 
modes as 



Y(x) = 
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(30) 

where we also imposed Y(x -C 1) ~ x to reproduce the 
correct thermodynamic limit. The chord length has only 
the first mode and so corresponds to kj = for any j. 
This expansion in terms of Fourier modes is particularly 
useful, because we expect that the contribution of the 
first few modes will be enough to have a reasonable ap- 
proximation of the scaling function Y(x). Indeed, Fig. [3] 
shows that only the first term k\ is enough to describe ac- 
curately the observed behavior for the RSP entanglement 
entropy 
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(31) 



with ki w 0.115. The obtained scaling function in pres- 
ence of disorder is greater than the chord length. 

Fig. [2] shows that the finite-size scaling in the ab-initio 
calculation are equivalent to the RSP ones (else for £ ~ L 
the various curves should bend). This means that the 
finite-size scaling of all S a in the spin chain is described 
by Eq. (31 1, as we also checked directly. 



C. Probability distribution of the Renyi entropy 

The disorder averaged Renyi entropy S a gives only ac- 
cess to the averaged number of the in-out singlets, while 
S a gives access to the full in-out singlets distribution 
P(n), i.e. the probability distribution of the Renyi en- 
tropy and so to the full entanglement spectrum. Indeed 
S a is related to the the cumulant generating function g(t) 
of the in-out singlets distribution by Eq. ( 16 1. 



We first consider the RSP data, because they allow to 
explore larger system sizes. Only after having established 
the asymptotic behavior we will consider ab-initio data 
and show consistency with the proposed scaling. 

We observed that the Renyi entropies 5^ SP do not 
have subleading corrections depending on the parity of 
the block, making the asymptotic analysis quite straight- 
forward. Oppositely, the data for S^ SP (see Fig. ^ show 
that they depend on the block parity in a way similar to 
clean systems^!. To analyze the numerical data we con- 
jecture the following asymptotic behavior 



S« SP (£) « 
A t S* SP ( 



B t \n2-(~l) e f t (S* SP (£))ln2, (32) 



where t is defined in Eq. ( |T7|. A t and B t are the two 
functions introduced in Eq. (|23|, while f t takes into ac- 
count the corrections to the scaling and goes to for 
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Figure 4: Top: RSP results for S a as a function of S a for a 
chain of 1024 spins and 1440000 disorder realizations. Bot- 
tom: Even-odd average of S a eliminating leading corrections 
to the scaling. In both panels, the continuous lines are the 
analytic RSRG result for A t . 



are linear only for a close to 1). 

From this linear dependence we can extract the func- 
tions At and B t using the RSRG relation 



(34) 



The resulting values for the universal coefficient A t ( a ) for 
a < 10 and for L = 1024 and L — 10000 are reported in 
Fig. [5] For small a (< 3.5) there are negligible fiinite- 
size corrections and the data perfectly agree with the 
RSRG result in Eq. (24 1, showing the predictive power 
of the RSRG to determine A t . For larger a, finite-size 
corrections are important and indeed the data differ from 
the analytical prediction, but the larger system sizes are 
closer. We believe that in the thermodynamic limit the 
RSRG A t describes the correct behavior for any a. The 
reason of these finite-size effects is also easily understood: 
the asymptotic formula is valid for S a large, while in this 
region of a we have S a ~ 1. Even if not asymptotic, 
the large a results show an interesting behavior: inde- 
pendently of L, they follow a — 1/t behavior (see inset in 
Fig. [5]), typical of a Poissonian distribution of singlets. 
The reason of this Poissonian behavior can be traced back 
to the fact that for t — > — 00 we are giving a large weight 
to short-range singlets that are produced almost inde- 
pendently. Little weight is instead given to long-range 
singlets responsible of the universal physics and so for 
these values of a and L we are probing the UV physics. 
According to this interpretation, a crossover from the uni- 



versal behavior of Eq. ( 24 ) to a UV Poissonian behavior 
always takes place for a ~ In L, in agreement with Fig. 

m 

We can now move to the ab-initio calculation to check 
the validity of the RSP scenario for S a . As before, we 



I — >• 00. The form of the corrections is inspired by the 
results in clean systems, while the leading term is the so- 
lution asymptotic g(t) in Eqs. (23 1 and (24 1. In the top 
of Fig. |4]we also report the RSRG value for A t that seems 
to be in qualitative agreement with the numerical data. 
A full quantitative description requires the elimination of 
the corrections to the scaling. 

In order to provide an unbiased description of the 
asymptotic behavior of S a , we define the functions 
s&v en (£) anc j g odd^ f rom the interpolation relative to 
even and odd blocks respectively. We can isolate the 
leading behavior of S PSP by considering the average over 
the two interpolating functions, i.e. 




(33) 



This definition eliminates the leading corrections to the 
scaling. In fact, in the lower panel of Fig. [4] we have 
a linear relation between and S BSP for all reported 
values of a (while the non-averaged data in the top panel 



Figure 5: The universal constant A t obtained from RSP data 
for L = 1024 (1440000 disorder realizations) and L = 10000 
(320000 realizations). Main plot: For a < 3.5 finite-size ef- 
fects are negligible and the RSRG prediction (continuous line) 
describes the data. Inset: Crossover to the non-universal Pois- 
sonian behavior (green continuous line) for larger a. 
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Figure 6: Ab-initio S a as a function of S« for a spin-chain 
of 1024 spins and 73000 disorder realizations. The continu- 
ous lines represent the RSRG prediction for the slope. The 
additive terms are different from those in Fig. [4] 



Figure 7: Sc alin g functions for the correction to the scaling 
ft (S) in Eq. Q obtained as difference of s° dd (£) and s° vcn {I) . 
Full and dashed lines correspond to uniform and exponential 
distributions of disorder respectively. 



focus on the relation between S a and S a and in particu- 
lar on the universal slope of the linear relation between 
them. The results are reported in Fig. [6] Asymptotically, 
the slopes of these curves tend to the RSRG prediction 
for At shown as continuous lines in the figure. Also the 
finite-size scaling scaling is well described by Eq. (311, 



as evident from the fact that the linear relation between 
S a and S a is correct even for large values of I (i.e. of 
S a ) in the various plots. However, as clear by a visual 
comparison between Figs. [6] and [4] (top) , the constant 
term in this relation is different (and both different from 
the analytic B t in Eq. (p4|). The degree of universality 
of this term is discussed in next subsection. 

Having established the correct asymptotic behavior we 
can consider the oscillating corrections to the scaling de- 
fined in Eq. ( 32 ) . The numerical estimate of ft (S) can 
be obtained as 



MS) 



^odd 



V) 



(35) 



where the dots denotes subsubleading terms (we recall 
^odd/even interpolations and so defined for any £). The 
data obtained in this way are reported in Fig. [7] The lin- 
ear behavior in log-scale shows that for a < 5 (for larger 
a further sub-leading corrections must be considerecP^l) 
ft(S a ) decays exponentially 



ft(x)=F t e- 



(36) 



i.e. a power-law correction in I. Vt( a ) is a new uni- 
versal critical exponent governing the corrections to the 
scaling of S a , analogous to the one introduced in clean 
systenuPEI We can see that v t ( a ) decreases with in- 
creasing a, but a precise numerical estima te is diffi- 
cult. For clean systems it has been showrP^Ml that 
v t(a) — 2K/a with K an a-independent exponent equal 
to the scaling dimension of a relevant operator. We can 



rule out this form for the random spin-chain, but the 
accuracy of our results does not allow to establish nu- 
merically an exact formula for the a dependence of the 
exponent. We also mention that the corrections to the 
scaling are of the same form also in ab-initio calcula- 
tions, as qualitatively clear from Fig. [6] and quantita- 
tively checked but not reported here. This shows the 
correctness of the RSP description and also that the real 
spin-chain does not introduce new leading corrections to 
the scaling in addition to the RSP ones. 



D. Universality 

All the results presented until now, both ab-initio and 
RSP have been obtained for random distributions of the 
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Figure 8: S a for two disorder distributions. The RSP data 
are for chains of 10000 spins and averaged over 320000 con- 
figurations. 
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Figure 9: S a as function of S a for a — 2.9 and for two disorder 
distributions (RSP data with L = 10000 and 320000 configu- 
rations). The scaling function is disorder independent. 



Figure 10: The quantity A defined in Eq. (38 1 vs S a for uni 



form and exponential distributions of disorder. With varying 
a the two differences are the same, showing the universality 
of the coefficient B t ( a y 



coupling constant J in the interval [0, 1]. However, the 
universal prediction of RSRG must be independent of 
the distributions of J (as long as new symmetries are 
not introduced). We check this universality by study- 
ing the RSP chain with L = 10000 spins with coupling 
distributed both uniformly J G [0,1] and exponentially 
P(J) ~ e~ J . In Fig. [8] we report the numerical RSP 
results of S a for a = 1, 2 for the two distributions of 
the disorder. As expected, the two distributions lead 
to slightly different results: only the leading logarithmic 
term in £ is universal, while the additive constant term 
is not. 

To check the universality of the leading term, Fig. [9] 
reports S a as function S a for a = 2.9 (other values of a 
leads to equivalent plots) for the two distributions. The 
two curves perfectly coincide, despite when they are plot- 
ted as function of I they are different. This means that 
all the non-universal behavior of the additive constants 
is washed out and we are left with a universal function. 
At first this result can seem surprising, but it is easy to 
realize that, in this kind of plots, the dependence on the 
non-universal cut-off or lattice scale a disappear and the 
leftover difference of non-universal additive constants is 



universal. For example, for the conformal entropies (14) 
we have the universal relation 



s: 



CFT 



*j ■ 



CFT 
VN 



1 



where evidently all the a dependence disappeared. To our 
knowledge, this property has not been explored at all in 
clean systems, but e.g. one can easily check that in the 
exact results for the critical XY modepS, the dependence 
on the irrelevant parameter 7 disappears in Eq. ( 37 1 . 



Having established that both A t and B t are universal, 
we reconsider our results for the disordered systems. We 
already discussed for the uniform distribution (see Fig. 
[5]) how the numerical value of A t agrees with the ana- 



lytical RSRG prediction. The independence of A t on the 
disorder distribution confirms its universality. In Fig. [10] 
we plot the quantity 



1 - a 



(38) 



where gr^ (t) is the function in Eq. ( 22 1 and fi is fixed by 
S a via /1 = j^Sa + 3. This quantity has been built in 
such a way to cancel the leading behavior A t so to leave 
only B t . Albeit little noisy, Fig. [10] shows clearly the 
disorder independence of B t . 

Disappointingly, as shown for uniform disorder, we 
found that the RSP and ab-initio calculation for S a pro- 
vide different values for the constant B t i a \ that are both 



different from the RSRG expression in Eq. (24 1. On one 



hand, this is showing that the RSP description is unable 
to catch this feature of the spin-chain because numeri- 
cal RSP and ab-initio data disagree. On the other hand, 
this is also showing that while carrying out the analytic 
results for g(t), some of the assumptions made influence 
significantly this quantity. There are two possible expla- 
nations to motivate the last discrepancy. One is that the 
distribution /(/Lt) in Eq. ( 18 1 contains some additional 
(subleading) terms not considered here. In fact, as al- 
ready discussed, Eq. (18 1 has been deduced neglecting 



terms coming from the starting disorder distribution and 
it is only asymptotically true. The other possibility is 
instead that the discarded terms in the renewal equa- 

Several pieces of information 



19 1 contribute to B t 



tion 

have been indeed ignored there: memory beyond first or- 
der, multiple decimations, the flow of the distribution to 
the critical point, etc. We found it rather improbable 
that f(n) should be modified. It is difficult to imagine 
how to modify it keeping all the other correct results (i.e. 
the entanglement entropy, A t etc.). On the other hand, 
solving the renewal equation in the presence of the dis- 
carded effects is very hard (maybe impossible). Thus, to 
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convince ourself that these processes can be responsible 
of a changing in B tl we tried to add some oversimpli- 
fied processes (but physically motivated) to the renewal 
equation: we found that all these processes change B t , 
but leave A t unchanged, showing that this is the most 
probable explanation of the discrepancy. However, from 
the ab-initio results, we know that the real spin-chain 
introduces further corrections to this term and so we do 
not find reasonable to embark in a difficult calculation, 
that in any case will not provide the correct answer for 
the spin-chain. 

To conclude the universality section, it is worth to men- 
tion that the oscillating corrections to the scaling (the 
function f t in Eq. ( [32] ) ) also do not depend on the dis- 
order distribution as shown in Fig. [7] confirming their 
universality 



V. CONCLUSIONS 

We provided an analytical and numerical description of 
the Renyi entropies S a and S a in a random singlet phase. 
For S a the leading logarithmic behavior is a-independent 
and only the subleading constant term depends on a: 



written in a completely universal form as 



1 



1 



-InTr 



In 2 



ln£ + d' 



(39) 



The leading universal term has been determined analyt- 
ically, while the non-universal correction d' a only numer- 
ically. Oppositely, the leading universal term of S a has 
a non-trivial a-dependence. Its scaling behavior can be 



S a = i InTr [p a ] ~ An_ a n n2 S a + Bn_ a \ ln2 , (40) 

1 — a 

Indeed, we pointed out that the functions A t and bt con- 
necting linearly S a and S a are both independent of the 
cut-off length introduced by the chain, and so univer- 
sal. The analytic result based on the solution of real- 
space renormalization group equations agrees perfectly 
with the numerical data as shown in Fig. [5] giving a 
full characterization of the asymptotic behavior. Instead 
a first-order RG prediction for the subleading term B t 
disagrees with the numerical data. We showed evidences 
that this disagreement should be related to the approx- 
imations done in the RG equations. Only an improved, 
but much more difficult (and maybe impossible) calcula- 
tion can provide the exact result for B t . 

We also studied the finite-size scaling: for finite chains 
the above relations still hold if the subsystem length £ is 
replaced by a modified chord length that is phcnomcno- 
logically well approximated by Eq. (31 1. We do not have 



a theoretical explanation for this finite-size scaling form. 

Assuming that the random-singlet description is 
equally valid for the random Heisenberg model, as is 
plausible and often assumed but not yet proved or firmly 
confirmed numerically, then we are in the surprising situ- 
ation of knowing the entanglement spectrum exactly for 
the random Heisenberg model but only approximately 
for the correspondin g p ure model (apart from some ex- 
act results for small .P^l). 
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